MOMENTS OF ASKEY- WILSON POLYNOMIALS 
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^v^j ■ Abstract. We give new formulas for the moment fi n (a, b, c, d; q) of Askey- Wilson polynomials. 

As a corollary we obtain a symmetric polynomial expressions for fi n (a, b, c, 0; q). We give a 
1 combinatorial proof of the formula for /i n (a,b,0,0; q). We also give the first combinatorial 

proof of the formula for the moments of q-Laguerre polynomials due to Corteel, Josuat- Verges, 
Prellberg, and Rubey. We show that our formula can be used to derive various results in the 
_-. ■ literature. Two positivity conjectures are given. 

I— > 

1. Introduction 

. The monic Askcy- Wilson polynomials P n = P n (x;a,b,c,d;q) are denned by the three-term 

[ recurrence P n +\ = (x— b n )P n — \ n P n -i with P_i = and Pq = 1 for b n = g(a + a -1 — (A n + C„)) 

and A„ = jA n ^iC n , where 

(1 - abq n )(l - acq n )(l - adq n )(l - abcdq™- 1 ) 

CZ' A n - 



a ( 1 - abcdq 2n ~ 1 ) ( 1 - abcdq 2n ) 
a(l - q n )(l - bcq 71 - 1 )^ - bdq n - r )(l - cdq^ 1 ) 



(1 - abcdq 2n - 2 ){\ - abcdq 2 "- 1 ) 

J> ' We refer to [8] for the standard basic hypergeometric notation and for information about the 

\& , Askey- Wilson polynomials. 

The n th moment /j, n (a,b,c,d;q) of the measure w(x;a,b,c,d;q) for the Askey- Wilson polyno- 
^5" . mials 

j. |, n n (a,b,c,d;q)= j x n w(x;a,b,c,d;q)dx 

is a rational function of a, b, c, d and q. (We will not need the explicit expression for w(x; a, b, c, d; q), 
see [HI (6.3.1), (7.5.15)].) This moment fi n (a, b, c, d; q) can be explicitly given as a double sum (see 
[1 Theorem 1.12]) 

^_ 1 >p {ab,ac,ad;q) m m ^ g- j2 q~ 2j (qg J + q~ 3 /a) n 

^" {) Vn{a,b,C,d ]q) - 2n ^ {ahcd . q)m 1 2^^ql-2 l/a 2. qUq . q 2 ] + l a 2. q)m _ 3 - 

<&'. 

However this expression is not obviously symmetric in a, 6, c, and d, even though the polynomials 
and moments are symmetric, nor does it exhibit the correct poles of fi n (a,b,c,d;q) as a rational 
function. 

The purpose of this paper is twofold. First we give new expressions for the moments \x n (a, b, c, d; q) , 
which are symmetric and polynomial when d = (see Theorem II .2j) . We also give a symmetric 
version for all a, b, c, d in Theorem 12. 5[ although the polynomial dependence in q is not clear. We 
give new expressions for the moments /x ra (a, 6, c, d; q) in the special case b — —a, d = — c, Theorem 
2.31 and Theorem II .61 We give two new positivity conjectures, Conjecture [T] and Conjecture O 

Our second goal is to combinatorially study the moments ^t n (a, b, c, d; g) as functions of a, b, c, 
and d. Corteel and Williams [7] give a combinatorial interpretation for the moments /x„(a, 6, c, d; q) 
using a rational transformation over the complex numbers of the parameters a, b, c, and d to 
parameters a, /?, 7, and 5. (See Section|4]). Using their ideas explicit coefficients of certain terms in 
Hn(a, b, c, d; q) as Catalan numbers are given in Theorems l4.2l and l4.6l The moments ^ n (a, b, c, d; q) 
are the generating function for Motzkin paths with weights which are rational functions of a, b, c, d 
and q. We use this setup and a generalization of an idea of D. Kim |16j to combinatorially prove 
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Theorem 11.31 and Corollary 11.41 in Section [3] Thus we give the first combinatorial proof of the 
formula of Cortccl et al. for the moments of q-Laguerre polynomials that is essentially equivalent 
Corollarv ll.4[ see Section[5] Josuat- Verges [13] gave a different combinatorial proof of Theorem ll.3[ 
but our proof is the first combinatorial proof of Corollarv ll.4l 

We now state the explicit formulas for the moments /i n (a, &, c, d; q). 

Using (fTJ) we prove the following theorem, which exhibits 2 n (abcd; q)„/i„(a, b, c, d; q) as a sym- 
metric polynomial in b, c, and d. 

Theorem 1.1. The Askey-Wilson moments are 



2 n fi n (a,b,c,d;q) = £ 



(ab, ac, ad; q), 
(abed; q) m 



n+l 



■h) m E 



s-l 



n — 2s — m 

x £ a 

p=0 



n+2s+2p 


m + p 




n 


- 2s - p 




m 


q 




m 



(-n+2s+p)m+(™) 



When d = 0, we obtain an explicit formula which is a polynomial in each parameter a, b, c, 
and q. 

Theorem 1.2. The Askey-Wilson moments for d — are 
2> n (a, b, c, 0; q) = ( («"fc) - x ) ) 

J2 a u b v c w (~l) t q^ 1 ) 



k=0 

X 



U + V + t 




V + w + t 




w + u + t 


V 


Q 


w 


Q 


u 



u+v+w+2t—k 

where the second sum is over all integers < u, v, w < k and — k < t < k/2 satisfying u + v + w + 
2t = k. 

The special case c = of Theorem 11.21 is equivalent to a result of Josuat- Verges [131 Theo- 
rem 6.1.1]. 

Theorem 1.3. The Askey-Wilson moments for c = d = are 
(2) 2> n (aA0Aq) = it((j^)-(nJ_ 1 )) £ a^-l)^ 



U + V + t 
U, V, t 



k=Q / / u +v+2t=k 

where the second sum is over all nonnegative integers it, v, t satisfying u + v + 2t = k. 

There are two interesting special cases of Theorem 11.31 which are the next two corollaries 
Corollary 1.4. We have 

k 



2 >„(a, q/a, 0, 0; ,) = £(( ^ ]-( ^_ J ] [q/af £ 1 



i=0 



Corollary 1.5. We have /i2n+i(a, — a, 0, 0; q) = and 



4™^2n(a,-a,0,0;q) = £ 



fc=0 



2n 
i — k 



2n 

n — k — 



J £(-ir? (T W) 



k-%a 



2k-2i 



i=0 



2k - i 
i 



lib = —a and d = — c, the Askey-Wilson measure is an even function, so the odd moments are 
0. In this case there are alternative expressions for the even moments, not obtained by specializing 
(fTJ) or Theorem O 
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Theorem 1.6. The non-zero Askey-Wilson moments for b - 

(-a 2 ;q) 2m (a 2 c 2 ;q 2 ) r 



-a and d 



-c are 



2n+2 



;<?) = 



rn— 



(ga 2 c 2 ; g 2 ) 



: (-<?T 



E 

s=0 



2n + 1 



2n + 1 
s - 1 



n — m — s 

E - 

p=0 



2n+4p+2s 


TO + p 




n — p — s 




m 


<? 2 


m 



-2m(n-p-s)+m(m-l) 



2. Askey- Wilson moments 



In this section we consider the moments ^ n {a, b, c, <i; q) as functions of the parameters a, &, c, d 
and q. Our goal is to give new explicit formulas for these moments, using simple series and integral 
evaluations. We shall prove the specific results in the introduction: Theorem 11.11 Theorem 11.21 
Theorem 12. 31 Corollarv ll.41 Corollary 1 1.51 and Theorem II .61 We use (J} and the method of proof 
of ([1]) which appears in [5] to prove these results. We generalize ([1} and Theorem II .11 to a version 
(Theorem [53]) which is explicitly symmetric in all of the parameters a, b, c, and d. Two positivity 
conjectures, Conjecture Q] and Conjecture [2J are also given here. 

First we note the polynomial behavior of the moments. An explicit combinatorial version of 
this proposition, although with complex weights, is given in Proposition |4T] 



Proposition 2.1. 2" (abcd; q) n /i n (a, b, c, d; q) is a polynomial in a, b, c, d, q with integer coefficients. 

Proof. By ([1]), 2 n (abcd;q) n p n (a,b,c,d;q) is a polynomial in b, c, and d with possibly rational 
function coefficients in a and q. By symmetry in a, b, c, and d, it is also a polynomial in a, thus 
a polynomial in a, b, c and d, with coefficients which are rational function of q. The poles of 
these rational functions must lie on either \q\ = 1 or q = 0. On the other hand p n (a, 6, c, d; q) is a 
polynomial in the three-term recurrence relation coefficients b n and A n of Section [H with integral 
coefficients. Since b n and X n have no poles at these locations, 2 n (abcd; q) n p n (a, 6, c, d\ q) is also a 
polynomial in q, with integral coefficients. □ 

We proceed to the proofs of the results given in the Introduction. 

Proof of Theorem IJ.il We shall show that |T]) implies Theorem 11.11 

Apply the binomial and g-binomial theorems to find that the j-sum of |lj equals 



3=0 
m 

: E< 

j=0 



(-l)V 



{a 2 qi;q) m+1 (q;q) r 



-iy q 



E 

1 s=0 



n 
s-1 



a 2 V s E 

p=0 



TO + p 
TO 



a 2 Pq 3P 



Thus the j-sum is summable by the g-binomial theorem to obtain 



2 n /z„(a, b, c,d; q) 



(a&, ac, ad; q) m ^ m ^ 



m=0 



(<7, a&cd; g) r 



s=0 



n 
s- 1 



p=0 



m + p 
rn 



-n+2s+p , 



We extend the p and s sums from — oo to oo by declaring these extended binomial and q- binomial 
coefficients to be zero. Replacing sbyn+l-s, and then p by p — n— 1 + 2s for the second term 
of the difference gives 



on i u a \ {ab 1 ac,ad;q) m ^ 



(q, a&cd; g) 



E 



-n+2s+2p 



TO + p 

m 



s=0 



(3 



— n+2s+p , 



q) m - 



p — n — 1 

TO 



2s 
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The first term is zero unless p > and (—71 + 2s + p < —to or — n + 2s + p > 1). The second term 
is zero unless p — n — 1 + 2s > and {p + 1 < — m or p + 1 > 1). If p > and — n + 2s + p > 1, 
then the first and the second terms are nonzero and equal. Thus the difference is zero unless 
< p < n — 2s — to for the first term and < p — n— l + 2s and p + 1 < —to for the second term. 
Undoing the change of summation variables for the second term then gives Theorem 11.11 □ 



Proof of Theorem \1.SX We show that Theorem 11.21 follows from Theorem 11.11 

In Theorem 11.11 expand the terms (ab; q) m and (ac; q) m as polynomial in a, b and c by the 
g-binomial theorem to obtain 



2 n n n (a,b,c,0;q) 



n+l 

E 

s=0 



where 



x= V(-i) 



u.v>0 



-l) u + v q(t)+(l)- 



n 
s- 1 



n— 2s— p 

• E 

m=0 



n-2s 

E« 

p=0 



-n+2s+2p 



(g p+1 ,g-™ +2s+p ;g) m g r ' 



(q;q)u(q;q)v ^ Q (q;q) m - u (q;q) r 
The TO-sum is summable by the g-Vandermonde identity [5J II. 6] and we get 



X = (-l) n - 2s -P- u - v q( T ' 



1 )a u+v b u c v 



2s 

v 



p + u 
u 



By replacing p and s with new summation variables t 
and k = n — 2s, we obtain Theorem 11.21 



■2s— p — u — v , w 



p + v 
n — 2s — p — u 



-n+2s + 2p+u+v, 
□ 



Theorem 11.31 follows immediately from Theorem 11.2 



The moment 2 n /j, n (a, b, 0, 0; q) in Theorem 11.31 is not a polynomial with positive coefficients. 
However, if we use a factor that is a power of 1 — q, we can express this moment using polynomials 
with positive coefficients. To accomplish this we need some definitions. 

Let Ai(n, m) denote the set of matchings on {1, 2, . . . , n) with to unmatched points. Note that 
A4(n,m) = unless n = m mod 2. For ir € A4(n,m) we define cr*(-7r) to be the number of 
pairs (ei, e-i) of elements ei, ei £ it such that e\ = e2 = {12, 32} with i\ < i-i < j\ < 32 or 

e i = e 2 = {32} with ii < i 2 < ji- Josuat- Verges [TH Proposition 5.1] showed the following 

(see also |H Proposition 15]): if n = to mod 2, we have 
(3) 



E * c 

ttG (n,m) 



1 



(1 - (7)(»-™)/ 2 ^ \ w-^- - fc 



- - 1 



(-l)^(^ 1 ) 



k + to 

fc 



For any polynomial /(x) in several variables x, let [x n ]/(x) denote the coefficient of x n in /(x) 
By ([2]) and ([3]), we get the following proposition. 



Proposition 2.2. We have 
2 n n n (a,b,0,0;q) 



u,v>0 



U TV 

a b 



u + v 
11 



(i-q) 



-v)/2 



E i 

1 TTGM(n,u+v) 

In particular, for integers n, u, v > with n = u + v mod 2 and n > u + v, we have 

[a u b v ]2 n fx n (a, b, 0, 0; q) = (1 - g)("-"-^/ 2 /(n, u + «; g), 

where f(n, u + v;q) is a polynomial in q with nonnegative integer coefficients. 

The moment /i„(a, b, c, 0; g) does not have the property given in Proposition 12.21 Indeed we 
have 

[abc}2 :i fi 3 (a, b, c, 0; q) = q 3 + 3q 2 + 3g - 1, 
[a6 2 c 2 ]2 5 /i 5 (a, 6, c, 0; g) = {q 2 + q + l){q 6 + 2q 5 + 4q 4 + 4q 3 + 4<? 2 - 2g - 3). 



In the next section we give a combinatorial proof of Theorem 11.31 Finding a combinatorial 
proof of Theorem 11.21 is still open. 
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Problem 1. Find a combinatorial proof of Theorem 11.21 

Corollary 11.41 and Corollary [L5] are special cases of Theorem II .3 



Proof of Corollary \H\ Put b = q/a in Theorem 11.31 and substitute v = k — 2t — u, and then 
t + u = i. The resulting i-sum is evaluable by the q-Vandermonde sum [8l II. 6]. □ 

Corollary 11.41 explicitly proves that 2 n fi n (a,q/a,0,0;q) is a Laurent polynomial in a, with 
coefficients that are polynomials in q, with positive integer coefficients. The following conjecture 
is its analogue for /j, n (a,q/a,c,q/c; q). 

Conjecture 1. For each non-negative integer n, [n + l] g !2 n /x n (a, q/a, c, q/c; q) is a Laurent poly- 
nomial in a, c, and a polynomial in q with non-negative integer coefficients. 

One can show that the sum of the coefficients in Conjecture [T] is 2 n (n + 1)!. 
Proof of Corollary \1.5[ We can rewrite ([2]) as 



E 

k=0 



n — k 
2 



n—k i 

2 1 



Lfc/2j 

E 

i=0 



2 n M„(a,-a,0,0;g) = 
Then we are done by the Gaussian formula Q] Theorem 10, p. 71 
(4) £(-1) 

i=0 



k — i 
k-2i 



k-2i 
J=0 



(-1)' 



k - 2i 
j . 



n 




i 









0, if n is odd, 

(?;9 2 )n/2, if n is even. 



□ 



Another special case of the Askey- Wilson polynomials has a different expression for the mo- 
ments. Consider b — —a and d = — c, so that the Askey- Wilson measure is symmetric about the 
y-axis. In this case b n = 0, so the odd moments are zero, and the 2n th moment has a shorter 
alternative expression. 



Theorem 2.3. The non-zero Askey-Wilson moments for b 



-a and d 



-c are 



^/^nK -a, c, -c; q) = 



(~a 2 ;q) 2m (a 2 c 2 ;q 2 ) 



m=0 



(qa 2 c 2 ; q 2 



a 4j g 2j (ag- 7 



(q 2 , a 4 g 2 +^ ; <7 2 ) m _, {q 2 , a" ! 9 2 ) 3 



Proof. We use the idea of the proof of [6l Proposition 3.1]. Let y = 2x 2 — 1 = cos(20). We need 
to find the coefficients A m in the expansion 



1 



— E 



a 2 e 2l °a 2 e- 2lt 



m=0 



This may be done using the q- Taylor theorem, see [TTJ Theorem 1.1]. The resulting Askey-Wilson 
integral is evaluated as in [BJ Lemma 3.2]. We do not give the details. □ 

Theorem 12.31 can be simplified in exactly the same way that (JXJ) implied Theorem II .11 to obtain 
Theorem 11.61 

There is a positivity conjecture for the moments fi2n{a, — a, c, — c; q). These moments are not 
polynomials, but Proposition 12.11 implies that 4™(a 2 c 2 ; g)2nM2n( a j — c, — c; q) is a polynomial. 
Half of the apparent poles of /i2n(a, — a, c, — c; q) do not occur. 

Proposition 2.4. TTie Askey-Wilson moments 

T2n{a 2 1 c 2 ) = A n (qa 2 c 2 ; q 2 ) n ^2n{a 1 -a, c, -c; 9)/ (1 - 9)™ 

are polynomials in a 2 , c 2 and q with integer coefficients. Moreover the sum of the coefficients in 
T 2n (a 2 ,c 2 ) is 2 2n {2n- l)(2n - 3) • • • 1. 



6 



JANG SOO KIM AND DENNIS STANTON 



Proof. After rescaling the Askey- Wilson polynomials, we have T2 n (a 2 , c 2 ) = (qa 2 c 2 ; q 2 ) n v 2m where 
V2n are the 2n th moments for the orthogonal polynomials defined by 



where 



Pn+l(x) = Xp n (x) - \ n p n -l{x), 

(1 + aV _1 )(l + c 2 q n - r ){l - a 2 c 2 q n - 2 ) . 

An = ~ o o — ox ~ o o — ~ \n\, 



(1 - a 2 c 2 g 2 "- 3 )(l - a 2 c 2 q 2n - 1 ) r " lq ' 
From Viennot's theory [231 . V2 n is the generating function for Dyck paths from (0, 0) to (2n, 0) 
with weights given by a product of \\s. So 

Q(a 2 ,c 2 ,q) 



V2n 



nr =1 (i 



for some polynomial Q with integer coefficients. By Theorem 12.31 v 2n has simple poles, as a 
function of c 2 , possibly only at (qa 2 c 2 ;q 2 ) n = 0. So (qa 2 c 2 ; q 2 ) n ^2n is a polynomial in a 2 , c 2 and 
q, with integer coefficients. 

For the final assertion, if q = 1, A n = rt(l + a 2 )(l + c 2 )/(l — a 2 c 2 ), which is the three term recur- 
rence coefficient for a rescaled Hermite polynomial. The 2n th moment for the Hermite polynomials 
is (2n- l)(2n-3)---l. □ 

Conjecture 2. The coefficients of T2„(o 2 , c 2 ) are non-negative integers. 

If q — 0, one may show that T2 n (a 2 , c 2 ) is a non-negative polynomial by a combinatorial method. 
The sum of the coefficients is 2 2 ™. It is a generating function for certain non-crossing complete 
matchings. 

Although simple, ([1} does not clearly demonstrate the symmetry or polynomiality of fJ> n (a, b, c, d; q) 
in all four parameters a, b, c and d. We next give such a formula, which generalizes Theorem ll.il 
Let A be an arbitrary parameter. Let 

8 W 7 (m) = 8 W 7 (A 2 /q; A/a, A/b, A/c, A/d, q-" 1 ; q; abcdq m ). 

Note that (aA, bA, cA, dA; q) m gW^im) is a symmetric polynomial in a, 6, c, d: 



(5) 



(aA, 6A, cA, dA; q) m s W 7 (m) = ^ 



(A 2 /q;g), 1 - A 2 ^ 1 



^(A 2 q™;q), 1 - A 2 /q 



x (a6cd) J (A/ a, A/b, A/c, A/d; q)j{Aaq> ,Abq> ,Acq> ,Adq j ;q) m - :j . 

Using Watson's transformation [8j (III. 17)] of an %Wj to a 4</>3, the following apparent rational 
function of A and q is in fact a polynomial in each of the parameters: a, b, c, d, A, and q. 



(6) 



(aA,bA, cA,dA;q) T 
{A 2 ;q) m 



iW 7 (m) = 

3=0 



(cd) J (A/c, A/d; q)j(a6; q), 



x (Aaq J , Abq j , cd; q) m -j ■ 

The next result gives a symmetric polynomial version for 2 n (abcd;q) n p, n (a,b,c,d;q). Theo- 
rem [23] is independent of A. 



Theorem 2.5. 



2 n fi n (a,b,c,d;q) = 



(aA, bA, cA, dA; q) m 
(A 2 , abed; q) m 
n— 2s— m 

x £ A- 

p=0 



n+l 



-q) m sW 7 (m)J2 



n 
s - 1 



n+2s+2p 


m + p 




n — 2s — p 




m 


q 


m 



m(-n+2s+p)+(™) 



Proof. We again use the q- Taylor theorem and follow the proof of [6l Proposition 3.1]. We expand 
x n in terms of the basis (Ae , Ae~ lB ; q)k- The resulting integral is a special case of the Nasrallah- 
Rahman integral [8] (6.3.9), p. 158]. An analogue of ([T]) is obtained. As before, the same steps 
with the binomial and q-binomial theorems yield the stated result. □ 
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If A = a, then Theorem 12.51 becomes Theorem 11.11 Theorem 12.51 has one defect: not all of 
the powers of q are positive due to the q m (- n + 2s +p) term. The individual terms are Laurent 
polynomials in q. 

If A 2 = g, the p-sum in Theorem l2.5l is evaluable by the g-Vandermonde sum [8j II. 6]. 
Corollary 2.6. If A 2 = q, 



n+l 



X A 



-n+2s 



s=0 

n + to + 1 — 2s 
2to+ 1 



n 
s - 1 



An explicit polynomial version of Theorem 11.21 which is unwieldy, may be given by expanding 
each of the eight factors in ([5]) by the g-binomial theorem. The resulting formula would appear to 
be too complex to be useful. 



3. Combinatorial proofs 

In this section we give combinatorial proofs of Theorem 11.31 and Corollary 11.41 
The main idea is as follows. We first interpret the moment 2 ra /Lt„(a, b, 0, 0; q) as a weighted 
sum of Motzkin paths. Then using Penaud's decomposition [50] we can decompose a weighted 
Motzkin path into a pair of paths: a Dyck prefix and another weighted Moztkin path. We map 
the new weighted Moztkin path to a new object: doubly striped skew shapes. These objects 
are a generalization of striped skew shapes introduced by D. Kim [16] in order to prove the 
moment formula for Al-Salam-Carlitz polynomials. We then find a sign-reversing involution on 
doubly striped skew shapes and show that the fixed points have a weighted sum equal to a q- 
trinomial coefficient. This completes the proof of Theorem 11.31 For Corollary 11.41 we find a 
further cancellation on the doubly striped skew shapes which leaves only one fixed point for given 
size. 

A Motzkin path is a lattice path inMxN from (0, 0) to (n, 0) consisting of up steps (1, 1), down 
steps (1,-1), and horizontal steps (1,0). We say that the level of a step is i if it is an up step or 
a down step between the lines y = i — 1 and y = i, or it is a horizontal step on the line y = i. 
A weighted Motzkin path is a Motzkin path in which each step has a certain weight. The weight 
wt(p) of a weighted Motzkin path p is the product of the weights of all steps. Note that the level 
of an up step or a down step is at least 1 and the level of a horizontal step may be 0. 

Let Mot ra (a, b) denote the set of weighted Motzkin paths of length n such that 

• the weight of an up step of level i is either q % or — 1 , 

• the weight of a down step of level i is either a&q 1-1 or —1, 

• the weight of a horizontal step of level i is either aq % or bq l . 

Then by Viennot's theory [23] we have 

(7) 2>„(aA0,0;<?)= E wt ( p )' 

P6Mot„(a,b) 

We define Mot* (a, b) to be the set of weighted Motzkin paths in Mot„(a, b) such that there is 
no peak of weight 1, that is, an up step of weight —1 immediately followed by a down step of 
weight —1. 

By the same idea that is a variation of Penaud's decomposition as in [T3] Proposition 5.1], we 
have 

w E *t( J, ) = E((A)-( E=4 n _ 1 )) E 

PGMot„(a,6) fc=0 ^ 2 ' ^ 2 ' ' PeMot* (a,f>) 

By ([7} and J5]), Theorem II .31 and Corollary II .41 can be restated as follows. 
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% — i 



1 \ 



l-L 

-1 




1 



Hi 




FIGURE 1. Converting a weighted Motzkin path to a doubly striped skew shape. 



Theorem 3.1. We have 



u + v + t 



]T wt(P)= £ o^t-i)*^ 1 ) 

PGMot*(a,h) u+D+2t=/s 

£ wt(P) = (g/a) fe ^a 2i g i(fe - i - 1) . 

PGMot*(a,ij/a) i=0 

We note that the second identity in Theorem 13.11 is equivalent to Proposition 5 in [5] which is 
used to prove the formula (ITBT) for the moments of g-Laguerre polynomials. Corteel et al.'s proof 
of [5] is combinatorial except for the proof of Proposition 5 in [S]. In this section we prove the 
above theorem combinatorially, thus providing the first combinatorial proof of (|16p . 

In order to give combinatorial proofs of the two formulas in Theorem 13.11 we introduce doubly 
striped skew shapes. These objects are a generalization of striped skew shapes introduced by 
D. Kim [16]. 

A doubly striped skew shape of size to x n is a quadruple (A, [A, W, B) of partitions fi C A C (n m ) 
and a set W of white stripes and a set B of black stripes with W n B = 0. Here, a white stripe is 
a diagonal set S of \j \x such that \/fi contains neither the cell to the left of the leftmost cells of 
S nor the cell below the rightmost cell of S, where a diagonal set means a set of cells in row r + i 
and column s + i for i = 1, 2, . . . ,p for some integers r,s,p. Similarly, a black stripe is a diagonal 
set S of A/fi such that A/fi contains neither the cell above the leftmost cell of S and the cell to 
the right of the rightmost cell S. We will call a cell in a white stripe (resp. black stripe) a white 
dot (resp. black dot). 

Let DSS(to, n) denote the set of doubly striped skew shapes of size mxn. We define the weight 
of (A, (i, W, B) G DSS(m, n) to be 



(9) 



wt a<b (X,fi,W,B) = a m b n {-l)W+\ B W x M-^ 



Kl/ab) 



w\ 



where ||W|| and ||S|| are the total numbers of white dots and black dots respectively. 
There is a simple correspondence between Mot£(a, b) and doubly striped skew shapes. 



Lemma 3.2. We have 



E wt ( p ) 

PeMot*(o,t) 



i=0 SSDSS(i,fc-i) 



22 wt ,6(5). 
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FIGURE 2. An example of the bijection p. The up steps of weight y and the down 
steps of weight 1 are colored red. 

Proof. It is enough to show that there is a weight-preserving bijection p : Mot£(a, b) — > Ui=o DSS(i, k— 
i). Let P G Mot£(a, b). We will construct an upper path and a lower path which determine two 
partitions A and p respectively. The upper path and the lower path start at the origin. For each 
step of P, we add one step to the two lattice paths as follows. If the step is an up step of weight 
abq 1 ^ 1 , we add a north step of weight a to the upper path and an east step of weight b to the 
lower path. If the step is an up step of weight —1, we add a north step of weight a to the upper 
path and an east step of weight b to the lower path, and we make a white stripe between these 
two steps, see Figure [TJ Similarly we add one step to the upper and lower paths for other type of 
step in P as shown in Figure [TJ Then we define p(P) to be the resulting diagram. See Figure [2] 
for an example of p. It is easy to see that p is a weight-preserving bijection. □ 

By Lemma 13.21 Theorem 13.11 is equivalent to the following proposition. 
Proposition 3.3. For k > 0, we have 

k 

(10) Yl E w w(s) 

t=o seDss(t,fc— o 

k 

(11) E W Wa(S) 
»=0 S6DSS(i,fc-i) 

In order to prove Proposition ^ . 31 we need D. Kim's sign-reversing involution. We now recall his 
involution. By a sign-reversing involution on DSS(m, n) we mean an involution / : DSS(m, n) — > 
DSS(m, n) such that if 5 € DSS(m, n) is not a fixed point, i.e. f(S) ^ S, then wt(/(5)) = - wt(S). 

Let DSS(m, n;a, /3, 7, S) denote the set of (A,/i, W, B) e DSS(m,n) such that \ = a, p — (3, 
W = 7, and B = 8, where a can be "— " and in this case there is no restriction on A, and it is 
similar for j3, 7, and 5. For example, DSS(m, n; A, — , 0, — ) is the set of (A,/x, 0,-B) € DSS(m, n) 
for all possible p and B. 

The following is proved by D. Kim [161 Theorem 3.2]. Since we will use his proof we include it 
here. 

Lemma 3.4. For a partition A C {n m ), we have 

Wt„, 6 (5)= J2 wt a , b (S) = a m b n q^. 

SeDSS(m,fi;A,-,0,-) SeDSS(m,n;A, 0,0,0) 

Proof. We construct a sign-reversing involution <f> : DSS(m, n; A, — , 0, — ) — > DSS(m, n; A, — , 0, — ) 
as follows. 

Let (A, /i, 0,-B) € DSS(m, n; A, — , 0, — ). Suppose p has r nonempty rows and the lowest row 
containing a black dot is row s. If there is no black stripe, then we let s = 0. If r = s = 0, 
equivalently p = B = 0, we define (A, p, 0, B) to be a fixed point. Otherwise, we define 0(A, p, 0, B) 
as follows. 

Case 1: If s > r, move the dots in the stripe containing a dot in row s (there must be a unique 
such stripe) all the way to the left and fill the cells containing these dots. 

Case 2: If s < r, we remove the leftmost cells of the last t rows in p and build a black stripe of 
size t ending in row r for the unique integer t for which this operation is possible. 



J2 a u b v {-iyq( t+ *) 

u+v+2t=k 
k 

(q/a) k J2a 2i 1 i{k ~^ 1} - 



u + v + t 

U, V, t 



i=0 
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Figure 3. Two examples of the involution 
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Figure 4. Slide rules 



See Figure [3] for some examples of this map. It is easy to check that <j> is a sign-reversing 
involution on DSS(m, n; A, — , 0, — ) with only one fixed point (A, 0, 0, 0), which implies the desired 
identity. □ 

The involution <f> of D. Kim cannot be directly applied to doubly striped skew shapes with white 
stripes. However, in the proof of the next lemma, we show that there is a simple way to extend 
his involution when we have white stripes. 



Lemma 3.5. For a partition A C {n m ), we have 

E wt a , 6 (5) = 

SeDSS(m,n;A,-,-,-) SeDSS(m,n;A,l 



Wt ,6(5). 



Proof. Let (A, /i, W, B) e DSS(to, n). Ignoring the white stripes, we apply the involution <fi. Then 
the newly added cells may overlap with white stripes. In this case we slide each overlapped white 
stripe one step to the right. If the rightmost dot in the stripe cannot be moved, then we delete 
this dot and the cell containing it from A and slide the white strip. See Figure 2] We define 
ip(\, fi, W, B) to be the resulting doubly striped skew shape. 

One can easily see that ip is a sign- reversing involution on DSS(m,n) with fixed point set 
DSS(m,n;A,0,— ,0). Thus we get the desired identity. □ 

Now we are ready to prove Proposition 

Proof of (fTTTj) . By Lemma 13.51 we have 



x\-\\w\\ 



i=0 SeDSS(i,fc-i) 



i=0 SeDSS(i,fc-i;-,0,-,0) 

= J2 <m-i)Y E i l 

u+v+2t=k (\,<!>,W$)£DSS{u+t,v+t),\W\=t 

It remains to show that for fixed nonnegative integers u,v,t we have 

(\,@,W,®)£DSS(u+t,v+t),\W\=t 



(12) 



U, V, t 
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Figure 5. An example of the involution ip. 
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Figure 6. An example of the intermediate step of the map from DSS (it + i, v + t) 
to 5(0", 1', 2 V ). The word corresponding to this example is 2021000202211. 



Suppose (A, 0, W, 0) € DSS(u + t, v + t) and \W\ — t. For each white stripe we mark the row 
containing the leftmost white dot and chop off the column containing the rightmost white dot, 
see Figure [H Then the resulting partition is contained in a (u + 1) x v rectangle. Consider this 
partition as a lattice path from (0, 0) to (v, u + t) with three steps: an east step, a north step, and 
a marked north step. Here a marked north step corresponds to a marked row of the partition. 
Then we define the word w = w\W2 ■ ■ ■ w u + v +t by Wi = 0, 1, or 2 if the ith step is a north step, 
a marked north step, or an east step respectively. Then w is a permutation of u 0's, t l's and 
v 2's. It is not hard to see that (A, 0, W, 0) i-> w is a bijection from DSS(« + 1, v + 1) to the set 
S(0 U , 1', 2 V ) of permutations of u 0's, t l's and v 2's such that 



|A| - ||W|| = inv(w) + 



where inv(w) is the number of pairs of integers such that i < j and Wi > Wj. Since 



w£S(O u , 1*, 2") 



U + V + t 
U, V, t 



we get H 



□ 



Proof of (fTTjl . For S € DSS(m, n), let S" be the doubly striped skew shape obtained by rotating 
S by an angle of 180° and exchanging blacks stripes and white stripes. The map S h- » S' is clearly 
a bijection. If 5 = (A, fj,, W, B) € DSS (m,n), then 5' = (A', //, W, fl') € DSS(m,n) satisfies 
A' = (n m ) - fj,, fjf = (n m ) - A, |W'| = \B\, \B'\ = \W\, \\W'\\ = \\B\\, and ||S'|| = ||W||. Thus 
wt a>6 (5) = wt aifc (5')(g/afe) |w/|_|s| - In particular, we have wt 0ig / (5) = wt 0>g / (S")< 
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Figure 7. A staircase tableau and the labeling of its empty cells. 



By Lemma 13.51 the map S M> S' and Lemma 13.41 we have 



E E wt a , g/Q (5) = ^ 

i=0 SeDSS(i,fc-i) 



which finishes the proof. 



E wt 0)9 / (5) 

4=0 SeDSS(i,fc-i;-,0,-,0) 

E E wt Q)<?/a (5') 

i=0 S'eDSS(i,fc-i;((fc-i) i ),-,0,-) 
fc k 

E E wt ,, /0 (5') = EXte/ )* - ^* -0 . 

i=0 S'GDSS(i,fc-i;((fc-i) i ), 0,0,0) 1=0 



□ 



4. Staircase tableaux 

In this section we review the combinatorial model, called staircase tableaux, for the moments 
of Askey- Wilson polynomials The staircase tableaux were first introduced in [7] and further 
studied in [6\ Using the staircase tableaux we shall find the coefficient of the highest term in 
2 n (abcd; g)„/i„(a, 6, c, d; q). 

A staircase tableau of size n is a filling of the Young diagram of the staircase partition (n, n — 
1, . . . , 1) with a, /3, 7, S such that 

(1) every diagonal cell is nonempty, 

(2) all cells above an a or 7 in the same column are empty, 

(3) all cells to the left of a p or 8 in the same row are empty. 

Here a diagonal cell is a cell in the ith row and (n + 1 — i)th column for some i € {1, 2, . . . , n}. 
We denote by T(n) the set of staircase tableaux of size n. 

Each empty cell s of T € T(n) is labeled uniquely as follows. Here, for brevity let right(s) be 
the first nonempty cell to the right of s in the same row, and below(s) the first nonempty cell 
below s in the same column. 

(1) If right(s) has a P, then s is labeled with u. 

(2) If right(s) has a 8, then s is labeled with q. 

(3) If right(s) has an a or 7, and below(s) has an a or S, then s is labeled with u. 

(4) If right(s) has an a or 7, and below(s) has an j3 or 7, then s is labeled with q. 

See Figure [7] for an example of a staircase tableau and the labeling of its empty cells. 

For T € T(n), we define b(T) to be the number of a's and <S's on the diagonal cells, and A(T), 
B(T), C(T), D(T), E(T) to be the number of a's, /3's, 7's, S's, empty cells labeled with q in T 
respectively. For example, if T is the staircase tableau in Figure [3 we have b(T) = 3, A(T) = 2, 
B(T) = 3, C(T) = 3, D(T) = 3, and J5(T) = 11. 

Let 

Z n (y;a,P, 7 ,S; q )= £ y b ( T ) a A{T) pB(T) ^C(T) fiD(T) gE(T) ^ 
TeTO) 
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Figure 8. A staircase tableau without a's and /3's and the corresponding alter- 
native tableau. 



Corteel et al. [5] showed that 

(1 - q) n 

(13) Mn(a, c, d; g) = — ^„(-l; a, /3, 7, <5; <?)> 

2"^" n?=oW-7&Z J ) 

where 

1-9 o 1 - 9 ac (l - <?) j - q) 

a = — — — , p = — — — , 7 — 



(l + a^)(l + ci) , (1 - bi)(l - di)' (l + ai)(l + a)' (1 - bi)(l - di) ' 

Since 

a/3 — -jSq 3 



(1 -qf{l -abcdqi) 



(l + ai)(l + ci)(l-bi)(l-di) 
we can rewrite (|13[) as the next proposition. 



Proposition 4.1. The Askey-Wilson moments satisfy 

T(abcd;q) n ^(a,b,c,d;q) = i~ n £ (-1)^(1 ?) A(t) +B (t) +C (t) +D (t)-„ ? K(t) 

TeT(n) 

x (ac^M)^ ((1 + at)(l + Ci ))- A(T) - C(T) ((1 - W)(l - ^))- S(T) ^ (T) . 



The highest degree term appearing in Proposition 14.11 is a n b n c n d n q^ . This term can be 
tained when A(T) + B(T) + C(T) + D(T) + E(T) = ("J 1 ) and A(T) = B(T) = 0, and the 
3mcient is 

i -ti/_ 1 )b(T)/_ 1 ^C(T)+D(T)-n/_ 1 Nri-Cf(r)/_ 1 \n-D(r) _ f»(_l)b(T)_ 



Thus 

(14) [a"6™c"d"g(2)j 2 n (abcd;q) n ^ n (a,b,c,d;q) = « n (-l) b(T) > 

TeCT(ra) 

where CT(n) is the set of T e T(ra) with A(T) = B(T) = and C(T) + D(T) + E(T) = 
We will see below that the tableaux in CT(n) are effectively the same as Catalan tableaux in [53] . 
For the remainder of this section we will prove the following theorem. 

Theorem 4.2. We have 

a n b n c n d n q(*) 2 n (abcd; q)„fi n (a, b, c, d; q) = Cat , 

where Cat(7i) = -^j-j ( 2 r ™) if n is a nonnegative integer, and Cat(n) = otherwise. 

In order to prove the above theorem we recall alternative tableaux, which are in simple bijection 
with permutation tableaux, see [19j [25] . 

An alternative tableau is a filling of a Young diagram with possibly empty rows and columns 
with up arrows and left arrows such that there is no arrow pointing to another arrow. The size of 
an alternative tableau is the sum of the number of rows and columns. A cell is called a free cell if 
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there is no arrow pointing to this cell. An alternative tableau without free cells is called a Catalan 
tableau. Viennot |26| showed that the number of Catalan tableaux of size n is be equal to Cat(n). 

Suppose T € T(n) satisfies A(T) = B(T) = 0. Then nonempty cells of T have only 7 and 5. We 
replace every 7 with an up arrow and every 5 with a left arrow. Then by definition every diagonal 
cell contains an arrow and there is no arrow pointing to another arrow. For each diagonal cell we 
remove the column (resp. row) containing it if it has an up arrow (resp. left arrow). Then the 
resulting object is an alternative tableau. If b(T) = k then the corresponding alternative tableau 
has n — k rows and k columns, see Figure[8] If moreover T satisfies C(T) + D(T) + E(T) = ("J 1 ), 
thus T € CT(n), then the corresponding alternative tableau has no free cells, i.e., it is a Catalan 
tableau. Thus we obtain the following lemma. 

Lemma 4.3. The number ofTG CT(n) with b(T) — k is equal to the number of Catalan tableaux 
of size n with n — k rows. 

The following lemma is proved by Burstein [3j Theorem 3.1] using the permutation tableaux 
description. 

Lemma 4.4. There is a bijection between the set of Catalan tableaux of size n with k rows and 
the set of noncrossing partitions of {1, 2, . . . , n + 1} with k blocks. 

Since the number of noncrossing partitions of {1, 2, . . . , n} with k blocks is the Narayana number 
Nara(n,fc) = ^(fc)(fe-i)' b ^ Lemmas l4~3l and l4~4l we obtain 

n 

(15) «"(-l) b(T) =« n 5I(-l)' £ Nara(n + l,n+l-fc). 

TeCT(n) k=0 

There is a formula for the alternating sum of Narayana numbers. 
Lemma 4.5. [2j Proposition 2.2 corrected] 



(-l) fe Nara(n,fc) = (-1) 2 * 1 Cat (^Tp^ 



fe=i 

Now we can prove Theorem 14.21 
Proof of Theorem \4-2[ By (14]) and |15]) we have 

n+l 

a n b n c n d n q(z)] 2 n (abcd; q) n ^ n {a, b, c, d; q) = i n ^(-l)' i+1_ ' £ Nara(n + 1, fc). 

fc=i 

Then we are done by Lemma 14.51 □ 

By analyzing the Catalan tableaux one can prove the following. We will omit the proofs. 
Theorem 4.6. We have 



a n - 1 b n c n d n q& 



2 n (abcd] q) n [i n (a, 6, c, d;q) = — Cat 



2 



a n-i b n-i c n d n q ( n 2 )^ 2 n (abcd; q) n p n (a, b, c, d; q) = Cat {^^\ - Cat . 

5. Connections with other results 

Recently the problem of finding formulas for moments of orthogonal polynomials has gained 
some attention. In this section we will rederive several moment formulas using Theorem 11.31 and 
Corollaries 11.41 and 11.51 We will use the following fact. 

Lemma 5.1. // /i„ is the n th moment of the orthogonal polynomial P n (x) defined by 

P n +i(x) = (x- b n )P n (x) - A„P„_i(a;), 
and /j,' n is the n th moment of the orthogonal polynomials P„(x) defined by 
P; i+1 (x) = {x- c(d + b n ))P' n {x) - c^KP^x), 



MOMENTS OF ASKEY-WILSON POLYNOMIALS 



15 



then 



5.1. The g-Laguerre polynomials. The g-Laguerre polynomials are denned by b n = [n] q + 
y[n + l] q and A n = ylft] 2 ,. Let v n denote the n moment of the g-Laguerre polynomials. We refer 
the reader to 5, Proposition 1] for many interesting interpretations of v n involving permutations, 
permutation tableaux, matrix ansatz, etc. 

Josuat- Verges [14] and Corteel et al. [5] showed that 



n n—k 



v HJ k=0j=0 

Using Lemma 15. II we have 



3 + k 



n \ I n 
j -l)\j + k + 



1 ))E(- 1 )V<z 4(fe+1 - 



i=0 



: E (V 1/2 + y-^T- m (") 2™M m (-<Z2/ 1/2 , -2/" 1/2 , 0, 0; q). 

m=0 ^ ' 



By using Corollarv ll.4l for 2 m ^ m (— qy 1 ^ 2 , —y~ x l 2 , 0, 0; q), the binomial theorem for (y 1 ^+y~ 1 ^ 2 ) n ~ m , 
and Vandermonde's theorem to sum the m-sum, we obtain (|16|) . 

Using the rank generating function for the totally nonnegative Grassmannian cells, Williams 
[2"Tl Corollary 5.3] showed another formula for v n : 



(17) 



"« = E» i E(- 1 ) J [ i -^« 



i(j-i) 



n 
.7-1 



i=0 j=0 

We now show that the two formulas (|16|) and (JT7J are equivalent under simple manipulation. 
In other words, we will prove the following identity: 



as) (l-ffrEEc- 1 )^--^ 



iU-i) 



i=0 j=0 



n n — k 

k=a j=o 



n\ I n 
3 J \3 + k 



n 



n \ I n 
j - 1 J \j + k + 



i i(k+l-i) _ 



By the binomial expansion for (1 — q) n [i — j] = (1 — q l 3 ) n , the left hand side of (|18|) is equal 



to 



EE(-i)v^E(l)(-^) fc ((;)« 



i=o j=o 



i-i 



n i n+1 



n 



ft \ / n 



i=0 j=0 fc=0 

The summand is anti-symmetric in k and j. Thus we can let k range from i + 1 to n + 1. Replacing 
k H> fc + j + 1 and i ^ i + j and changing the order of sums we get the right hand side of (fT8|) . 

5.2. Josuat- Verges's formula. Josuat- Verges 13J defined a partition function Z n (a,b,y,q) and 
showed that (1 — q) n Z n (a, b, y, q) is equal to the n th moment defined by b n = 1 + y + (a + by)q n 
and A„ = y(l — q l )(l — abq 1 ^ 1 ). He found the following formula for Z n (a, b, y, q) 



1 " L 2 J 

(19) Z n (a,b iyi q) = -——Y, E H/) l <? m 



k=0 i=0 
n-k-2 

X 

j=o 



k + i 
i 



— K — Zl K 

E/ ^ \ {j){k+2i+j) ~ (j-l) (k+2i+j+S) EZ 



r=0 



a r (yb) 



k—r 
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By rescaling we get 

n 

(1 - q) n Z n (a, b, y, q) = y n ' 2 £ (y 1 ' 2 + y^y 



m— 



2 m f i m (ay- 1 / 2 ,by 1 / 2 ,0,0;q). 



Thus we get in the same manner as |T 



" n—k / / \ / 

i—n o— n \\J / \J 



71 W 71 
j - 1 J U + fc + 1 



fc=0 j=0 

x ^ a M 6' u y , ' +t (-l) t g(^ 1 ) 

U+D+2t = fe 

It is easy to check that the above formula is equivalent to (fTT)| . 



u + u + 1 

u, v, t 



5.3. Josuat- Verges and Rubey's formula. Josuat- Verges and Rubey [TS] considered the 
moment R n (a,b,c,d;q) of the orthogonal polynomials defined by b n = d + (a + b)q n and \ n 
(1 — q n )(c — abq™^ 1 ) and showed that 



,l.h 



k + 1 + 1 



(20) R n ( a ,b,c,d- iq ) = J2^-Tj2 

fc=0 ^=o 



„n—k—l -i2l—n+k 



M£(a,b,c; q), 



where 



fc>0 



cb x qz,q 
aqz 



q,bz 



Using pop Josuat- Verges and Rubey |15j derive moment formulas for various orthogonal polyno- 
mials. 

By rescaling we have 



R n (a,b,c,d;q) = c n/2 (cfc -1 



in 



fi m {ac 



-V 2 : bc~ 1/2 A0;q). 



Thus we get 
(21) 

n n 

R n (a,b,c,d;q) = J2H c ^ dn ' m O " (*^-i)) E aW(-l)V^ 

fc— ra—k u^-v-\-2t—k 

It is easy to check using g-binomial theorems that (|2"Tj) is equivalent to (|20p. 

5.4. The (i, g)-Euler numbers. Kim [18] defined the (£, q)-Euler numbers E n (t,q) by 

1 



u + v + t 

U, V, t 



Y J E n {t,q)z n = 



n>0 



[ l ]q[ l ]t,qZ 

[2] 9 [2] t:9 z 



1 



and showed that 

E. 



1 

v HJ fc=0 



2n 

n — k 



2n 
n — k - 



where [n]t j9 = (1 — tq n )/(l — q) and 
fc 3 

(22) T fc (*, g )=^^(-iy +i t 2 V 



j=0 i=0 



k — j 



k — i 
j - i 



k-i-1 
j-i - 1. 



Note that we have (1 - q) 2n E n (t,q) = 2 2n u 2n {^/=qt, s/^qE, 0, 0; g). Thus by Corollary O 
implies 
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(23) 
E n {t,q) = 



(1-9) 



2n 



E 

fc=0 



2n 
n — k 



n 

n — k 



2k- i 

i 



We now show that the above two formulas for E n (t,q) are equivalent. It is sufficient to show 
that 

(24) {-l) k q k2 T k (t, g- 1 ) = £ 9© (?! g 2 )*-^ 7 % 

In (I22p if we replace j i— > fc — j , interchange the z-sum and the j-sum, and replace j i— > j + i, 
then the j-sum is summable by the g-binomial theorem. Thus the left hand side of (1241) is equal 
to 

Lfe/2J f _ ., Lfc/2J 

E.2i „2i 2 -i/ 2i+l 2\ K 2 | V )2i+l„2! 2 +./ 2i+3. „2\ 

t q \q i q h-a . + > * 9 (,<? , <? M— 2*— 1 

i=0 1 J 1 i=0 

Lfc/2J r n , n .n Lfc/2J 



fc- i- 1 

2 



i=0 



2fc- 2i 
2i 



I 9 13)9 jfe-2i-l 



9 i=0 



2/t - 2i - 1 
2i + l 



which is equal to the right hand side of (|2~4")) . 

Zeng [28 showed that 



(25) 



s„(i, (? ) = i-"EE(- 1 )' 



m— i— 



q 2m-2in+i -n-i[2 TO ] t , ? ![2i + 

[2i],!![2m - 2p nr=o,W 2fc + 2i + 2 ]* 2 ^ ' 



where [2m] t , 9 ! = n 2 =i [*k<? and [2i] 9 !! = 111=1 [ 2fc ]<?- 

It is straightforward to check that Zeng's formula is equivalent to Theorem 12.31 when a 
and c = 0. 



-qt 
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